Compound and scale mixture of vector 
and spherical matrix variate elliptical 

distributions 

Jose A. Diaz-Garcia * 

Department of Statistics and Computation 
25350 Buenavista, Saltillo, Coahuila, Mexico 
E-mail: jadiaz@uaaan.mx 

Ramon Gutierrez Jaimez 

Department of Statistics and O.R. 
University of Granada 
Granada 18071, Spain 
E-mail: rgjaimez@ugr.es 



Abstract 

Several matrix variate hypergeometric type distributions are derived. The com- 
pound distributions of left-spherical matrix variate elliptical distributions and in- 
verted hypergeometric type distributions with matrix arguments are then proposed. 
The scale mixture of left-spherical matrix variate elliptical distributions and uni- 
variate inverted hypergeometric type distributions is also derived as a particular 
case of the compound distribution approach. 

Introduction 

Four classes of matrix variate elliptical distributions have been defined and studied by 



Fang and Zhand (|l990h . That m x n random matrix variate X is said to have a matrix 



variate left-spherical distribution, the largest of the four classes class of matrix variate ellip- 
tical distributions, if it density function is given by 



c(m, n) 



| S |™/2|0|m/2 



h (S" X / 2 (X - m)'® _1 (X - /x)£- 1/2 



where ft is a real function, c(m,n) denotes the normalization constant, S is an m x m 
positive definite matrix, this fact being denoted as X > 0, is an n x n matrix, > 0, 
and fj, is an m X n matrix. 
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When implementing frequentist or Bayesian methods, one may be interested in s itua- 
tions where X has a density function of the following form, see lFang and Zhand {l990) and 
Fang and ii (jl999l ). 



Isi^U Mg-^-Mye^tx-M)). (i) 

This fact is denoted as X ~ £CS mxn (p,^,®,h). This condition is equivalent to consid- 
ering the function h as a symmetric function, i.e. g : g(AB) = g(BA) for any symmetric 
matrices A and B. This condition is equivalent to that in which h(A) depends on A only 
through its eigenvalues, in which case the function h(A.) can be expressed as /i(A(A)), where 
A(A) = diag(Ai, . . . , A m ) and Ai, . . . , A m are the eigenvalues of A. Two subclasses of matrix 
variate elliptical distributions are of particular interest: the vector and spherical matrix 
variate elliptical distributions. For these distributions, A(A) = tr(A) and A(A) represents 
any function of eigenvalues of A, respectively. Note that vector matrix variate elliptical 
distributions are a subclass of matrix variate spherical elliptical distributions. Many well- 
know distributions are examples of these subclasses; one such is the matrix variate normal 
distribution. Other variants include vector matrix variate elliptical distributions e.g. Pear- 
son type II, Pear s on ty pe VII, Kotz type , Bessel and Logistic, among many others, see 
Gupta and Varga (|l993l) . Yet other are matrix variate spherical elliptical distr ibutions, e.g. 



Pearso n Type II, Pearson type VII and Kotz type among many others, see iFang and Li 
(1999) 



In the vectorial case, Muirheadl ( 1982 . p. 33) proposed a means of generating a family 



of vector variate elliptical distributions from a normal distribution. In th e general case, this 
i dea h as been extended to the matrix variate elliptical distributions by iGupta and Vargal 



(|1993l . pp. 78-79 and Section 4.1). The s ituation in which the spe cific elliptical distribution 
is a matrix variate normal is studied in IGupta and Vargal (|l993l . Chapter 4). Generically, 
distributions obtained by this procedure are termed scale mixture matrix variate normal or 
elliptical distributions. 



ip 



Arslanl (|2005l ) proposed the scale mixture of th e vector variate Kotz type distribution, also 



termed the t-type or generalised t distribution. Diaz-Garcia and Gutierrez-Jaimezl ( 20091) 



extend this idea to matrix variate vector and spherical Kotz type distributions using two 
approaches: scale mixture and compound matrix variate distributions. 

These forms of obtaining vector or spherical matrix variate elliptical distributions are of 



parti cular interest from a Bayesian point of view liJammalamad aka et aLl. 119871 ; iFang and Lit 
19991 ) and in the context of shape theory, see ICaro-Lopera et al. I (|2008l ). 



This paper introduces several families of matrix variate elliptical distributions. Section[2] 
gives some results on integration, using zonal polynomials. In terms of these results, various 
matrix variate hypergeometric type distributions are proposed, as particular cases. These 
include well-known distributions such as central and nonccntral matrix variate inverted 
gamma (inverted Wishart) distributions and matrix variate central and noncentral beta 
type II distributions. In Section [3l assuming a hypergeometric type distribution for the 
matrix parameter in matrix variate normal and matricvariate T distributions, several families 
of matrix variate elliptical distributions are obtained using the compound matrix variate 
approach. Section [6] introduces the scale mixture of a matrix variate elliptical distribution, 
which is derived as a particular case of the compound matrix variate approach, from where, 
all the results given in Section [3] can be particularised to this case. 
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2 Preliminary results 



Consider the following notation and definitions: the hypergeometric functions p F q with 
matrix arguments are defined by 



p F q (a!, . . . ,a p ;h, . . . ,b q ;Y) = ^^2 



k=0 K 



(Qi)k ■ ■ ■ {a p ) K C K (Y) 
(h) K ■ ■ ■ (b q ) K M ' 



where Y is a complex symmetric m x m matrix, C K (Y) is the zonal polynomial of Y of 
degree k, k = (fci, . . . k m ), k\ > ■ ■ ■ > k m > 0, k± + ■ ■ ■ + k m = k and oi, . . . , a p , b\ 
are real or complex constants, 



, b q 



(o)« = lj(o-(*-l)/2)* <> 



i=l 



with (x) n = x(x + 1) • • • (x + n — 1), (ir)o = 1- 
The multivariate gamma function is defined as 



l)/2], 



and 



T m [a, k]=tt 



_ m(m-l)/4 



ijr[a+ti-(t-l)/2|, 



whereT m [a,K] = (o)«r m [a] with Re (a) > (m-l)/2. see lKhatril (fl96l) and lMuirheadl (|l982h . 
Similarly, 



T m [a, -k] 



7T 



m(m — 1) /4 



[^rla-fe, - (m-*)/2] 



l(m— 



1)/4 nr[a-fc m+ i-i-(*-l)/2], 



where Re(a) > (m — l)/2 + k%. From 



we obtain that 



(-x) q = (-l)i(x-q + l) q = 
T m [a, -k] = 



(-iyr[x + i] 

r[a-g + l] ' 

(-i) fe r m [ a ] 



(-a + (m+l)/2) K 
The multivariate beta function is defined as 

P m [a,b} = f |Y| a -( m+1 )/ 2 |I m - Y| b -< m+1 )/ 2 (dY) = 

Jo<Y<I m 



r m [a]r m [b] 

T m [a + b] 



where Re (a) > (m - l)/2 and Re(6) > (m - l)/2, see|Herj (|l955i p. 480). 

Many distributions in multivariate analysis can be express ed in a form involving hyper- 
geome tric functions with matrix arguments, as considered by iRouxl (Il975l) . see also Press 



(1982, Section 6.6.3, pp. 170-171). These distributions contain as particular cases the cen- 
tral and noncentral gamma (Wishart) and matrix variate beta type I and II distributions 
and are termed matrix variate hypergeometric gamma (Wishart) type and matrix variate 
hypergeometric beta type I and II distributions. In particular, for matrix variate hyperge- 
ometr i c beta type II distributions, we obtain an alternat i ve ex pression to the one given by 
IRouxl (|l975h . based on the following lemma from iKhatril (|l966h . 
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Lemma 2.1. //R is any arbitrary complex symmetric m x m matrix, then 

I |Yr-(™ +1 '/ 2 |I m + Y|-(^)C K (YR) = Mi^M g<t( _ R) 
Jy>o (-6+ (to + 1)/2) K 

w/iere Re(a) > (to - l)/2 and Re(6) > (m - l)/2 + fci. 

Proo/. As given in lKhatril (|l966T ). □ 



If Y > has a matrix variate hypergeometric beta type II distribution, the argument 
of the hypergeometric function involved, in Roux's version, is (I m + Y) , whereas in the 
version based on Khatri's lemma, its argument is Y, as we see below. The importance of 
this fact is made apparent in the next section. 

The next result is obtained immediately from Lemma 12.11 



Corollary 2.1. Let R be any arbitrary complex symmetric to x to matrix, then 

I |Y| a -( m+1 )/ 2 |I m + Y\-^+ b \F q ( ai , . . . , a p ; h, . . . , b q ; RY) 
Jy>o 

= P m [a,b]p+iF q+ i(ai, . . . ,a p ,a;6i, ...,b q> -b+ (m + l)/2; -R), 
where Re(a) > (to - l)/2 and Re(b) > (m - l)/2 + k\. 

Proof. The final results are obtained using the fact that C K (6Y) = b k C K (Y), for a 
constant b. □ 

As a consequence of Corollary 12.11 we have the following alternative definition of the 
matrix variate hypergeometric beta type II distributions. 

Definition 2.1. Let H be any arbitrary complex symmetric m x to matrix. Then Y has a 
matrix variate hypergeometric beta type II distribution if its density function is 

/ Y (Y) oc |Y| Q - ( m+1 >/ 2 |I m + Y\-^ +b \F q ( ai , a p ; , b q ; HY), Y > 0, 

where the constant of proportionality is 

1 



b] p+ iF q+ i(ai, ...,a p ,a;bi,..., b q , -b+(m + l)/2; — B) ' 

with Re(a) > (to - l)/2 and Re(6) > (to - l)/2 + fci. 

Let Y be a positi ye definite to x to matrix and let us define P = Y _1 then (dY) 
|p|-(m+i)(dp)_ From lRouxl (|l975l . eq. (4.1)), we have the following result. 



Lemma 2.2. Let T and H be complex symmetric m x to matrices with Re(H) > 0. And 
assume thatY has a hypergeometric matrix variate gamma type distribution. Then P = Y _1 
has a matrix variate inverted hypergeometric gamma type distribution with the 
following density function: 

/ P (P) cx ctr{-EP- 1 }|Pr-( m+1 )/ 2 p F 9 ( ai , a p ; , b q ; TP" 1 ), P > 0, 

where the constant of proportionality is 



r m [a] p+ iFq(ai, . . . ,a p ,a;h, . . . ,b q ;T S 
and Re(a) > (to - l)/2. 
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Observe that if in Lemma 12.21 St = 0, then P has a matrix variate central inverted 
gamma distribution. And if p = 0, q = 1 and a = b\ 1 then P has a matrix variate noncentral 
inverted Gamma distribution. 

Similarly, from Definition 12.11 we have 

Lemma 2.3. Let 3 be any arbitrary complex symmetric mxm matrix. Then P = Y -1 has 
a matrix variate inverted hypergeometric beta type II distribution and its density 
function is 

/ P (P) cx |P| b -(™+D/2| Im + p\-( a +b) pFq ( ai> hq - HP- 1 ), P > 0, 

where the constant of proportionality is 

1 

P m [ a , b]p+iF g +i(ai, ...,ap,a;bi,..., b q , -b+(m + l)/2; — S) ' 
and Re(a) > (m - l)/2 and Re(6) > (to - l)/2 + fc x . 



The distribution in Lemma 12.31 contains as particular cases the matrix variate central 
an d noncentral inverted beta t ype I I distributions. 

Van der Morwe and Rouxl ( 1974 ) studied other families of matrix variate hypergeometric 



distributions, based on the 2-F\(a, b\ c; — Y) hypergeometric function with a matrix argument. 
From this, and taking the limits of parameters a, b or c, they obtain the central and non- 
central matrix variate gamma distributions, the cent ral matrix variate be t a dist ribution and 
the matrix variate normal distribution. As we see in lMathai and Saxenal (|1966l ) many other 
well-knows distributions can be obtained as particular cases of this distribution. Next, we 
introduce this distribution using the multivariate Mellin transform, Mathai ( 1997t h 

Let g(Y) be a function of the positive definite mxm matrix Y. The Mellin transform 
of g(Y) is defined as 

M g (Y) = f |Yr-(™ +1 )/ 2 3 (Y)(dY), 
Jy>o 

where Re(a) > (m — l)/2. 

Lemma 2.4. The Mellin transform of g(Y) — 2-Fi(a, b;c; — Y) is given by 

( |Yr-(" l+1 )/ 2 2J Fi(a,6; C ;-Y)(dY) = PrA^b - a]P m [a - a,c- a] ^ 
Jy>o P m i a ,c-a\ 

whereRe(a) > (m-l)/2, Re(a-a) > (m-l)/2, Re(6-a) > (m-l)/2, Re(c-a) > (m-l)/2 
and Re(c — a) > (to — l)/2. 

Pr oof. From the integral representation of 2-F1, see Herz ( 19551 eq. (2.12)) and Muirhead 
(|1982| . Theorem 7.4.2), 



Y) = [ lY^+V/^F^a^q-YXdY) 

JY>0 



r > 

1 



o-(m+l)/2 / ipia-(m+l)/2 



^m[ a i c — a] Jy>0 J0<R<I„ 

x \I m - R| c - Q -(" 1+1 )/ 2 |I m - YR|- fo (dR)(rfY) 

| R |a-(m+l)/2, I _ R , c -a-(m+l)/2 



@ m [ a i c — a] 7o<r<i„ 

x [ |Y| Q -( m+1 )/2|I m - YR|- fe (dY)(dR). 
Jy>o 
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Note that, |I m - YR| = |I m - R 1 / 2 YR 1 / 2 | = \I m - W|, where R 1/2 is the positive definite 
square root of R, such that R = (R 1 / 2 ) 2 |Muirheadl . [l98l Theorem A9.3, p. 588) and |Y| = 
|R- 1 / 2 WR- 1 /2| = |w||R|-\ with W = R^YR 1 / 2 . Then, (dY) = |R|-( m+1 )/ 2 (dW), 
from where 

M 9 (Y) 



fl m [ a > c — a] 7o<r<i, 

X 

f3 m [a,b-a] 



ftm[ a i c ~ a ] Jo<R<I 

f3 m [a, b - a](3 m [a - a, c - a] 



| R |a-a— (m+l)/2|j _ J^jc— a-(rn+l)/2 

f |Y| Q -( m+1 )/ 2 |I m - W\- b (dW){dY). 
Av>o 

| R |a-a-( m +l)/2|j m _ R|c-a-(m+l)/2^ Y ) 



□ 



P m [a,c-a] 

This was proved bv lVan der Morwe and Roux (1974), who took the limit when a tends to 
infinity in Lemma l2~4l and found the Mellin transform of the function g(Y) — iFi(b; c; — Y). 
Alternatively, the results can be obtained directly by integration, as we shown below. 

Lemma 2.5. The Mellin transform of g(Y) = iFi(b;c; — Y) is given by 

f lYr ( m+ i)/2 lFlib . c ._ Y){dY ) = ^MT^\c}T m [b-a} ^ 
J Y >o T m [6jr m [c-aJ 

where Re(a) > (m - l)/2, Re(6 - a) > (to - l)/2 and Re(c - a) > (to - l)/2. 

Proof. Noting that from IMuirheadl (|l982l Theorem 7.4.2, p. 264) 

1 



2 Fi(c- b,a;c;I„ 



ftm[ a i c — a] 7o<r<i„ 
[3 m [b-a,c-b] 



| R | C -b-( m +l)/2|j m _ R|&-«-(™+l)/2( dR ) 



And from the Kummer relation discussed in IMuirheadl (|1982L Theorem 7.4.3, p. 265 and 
Theorem 7.3.4), we have 

M g (Y) = [ lYp-^+^iFifccj-YXdY) 
Jy>o 

= f |Y|a-(»™+i)/2 i etr(-Y)iFi(c — 6; c; Y)(dY) 
iy>o 



T m [a] 2 F l (c - b, a; c; I m ) 

r m [a]r m [c]r m [6 - a] 



□ 



r m [6]r m [c - a] 

Now, from Lemmas El and ES] taking Y = 3 1/2 Y3 1/2 with (<fY) = |H|( m+1 )/ 2 (dY), 
we obtain the following. 

Definition 2.2. Let H be any arbitrary complex symmetric to x to matrix. Y is said to 
have a matrix variate generalised hypergeometric distribution if, 



1. Its density function is 

l H | Q /3 m [a,c- a] 



/ Y (Y) = 



P m [a,b - a\(3 m [ 



a — a, c — a\ 



|Y| a -( m + 1V2 2 Fi(a, b; c; -SY), Y > 0, 



with Re(a) > (to - l)/2, Re(a-a) > (m - l)/2, Re(b - a) > (to - l)/2, Re(c-a) > 
(to - l)/2 and Re(c - a) > (to - l)/2. 



G 



2. Or 

/ Y (Y) = ' B 'Tr m[ f!r m[ r " ] 1 \^r im+1)/2 ^(b; c; -BY), Y > 0, 
with Re(a) > (m - l)/2, Re(6 - a) > (to - l)/2 and Re(c - a) > (to - l)/2. 

Lemma 2.6. Let H &e any arbitrary complex symmetric to x m matrix. It is said that 
P = Y _1 /ias a matrix variate inverted generalised hypergeometric distribution if, 

1. Its density function is 

/ P (P) = |Bm>,c- a ] (^/^(a^-SP- 1 ), P > 0, 

Praia, 6 - aj/3 m [a - a, c - aj 

MftRe(a) > (to-1)/2, Re(a - a) > (m-l)/2, Re(6 - a) > (to-1)/2, Re(c-a) > 
(to - l)/2 and Re(c - a) > (m - l)/2. 

£. Or 

/p(P) = l B | a ^W^[^ -«] |p|-^( m+ l)/2 if1i(5;c; .gp-lj p > 

with Re(a) > (to - l)/2 ; Re(6 - a) > (to - l)/2 and Re(c - a) > (to - l)/2. 
Proo/. Follows from Definition EH taking P = Y" 1 with (dY) = |p|-("^+i) (rfp). □ 



3 Compound elliptical distribution of a random matrix 



In this section we propose several families of elliptical distributions based on an exte nsion 
to the matrix variate case of the vector case idea, introduced by iMuir headj dl98 2j), the 
approach Known as compound d istribution. This approach was used by IRouxi (|19Tlh and 
Van der Morwe and Rouxl (|1974[ ) for the distribution of a positive definite random matrix. 
In general, assume that the conditional distribution of 



X|P ~ ££S mXn (n, S 1 / 2 *^ 1 / 2 , 0, h), 



(2) 



with * : ^jn( m +i)/2 _> s Rm (m+i)/2 ; ^(p) > Q . where p > o has the distribution function 
G(P). Then X has a left-spherical elliptical distribution (compound distribution) with a 
density function given by 



;(m, n) 



| S |n/2|0|m/2 



h (*(P)- 1 S- 1 / 2 (X - A t)'0- 1 (X - /x)S" 1/2 ) dG(P) 
p>o l*(P)l" /2 



(3) 



where ^(P) -1 denotes the inverse of the matrix ^(P) (not the inverted function of >&(•)). 

Let us now consider two particular matrix variate left-spheric al ellipt i cal di s tributions, 
the m atrix variate norm al and the mat ricvariate T distributions, see lDickevi (|l967luBox and Tiao 
(|l972l . pp. 441-448) and lPressI (Il982l . pp. 138-141). 



4 Compound matrix variate normal distribution 

Recall that if X ~ N m xn{lJ', S, 0), its density function is given by 

' etr{-iE- 1 (X-/x)'0- 1 (X-/x)} 



(27r)"m/2|£|»/2|0|W2 
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Theorem 4.1. Assume that X|P has a matrix variate normal distribution, 

where P has a matrix variate inverted hypergeometric gamma type distribution. By Lemma 
its density function is 



g p (P) cx ctr{-HP- 1 }|Pr a -( m+1 )/ 2 p F g (a 1 , . . . , a p ; b 1} ... , b q] YP^ 1 ), P > 0, 

where the constant of proportionality is 



r m [a] p +iF g (ai, . . . , a p , a; &i, . . . , b q ; Y H 1 ) 



and Re(a) > (m — l)/2. TTien X /ia,s a matrix variate left- spherical elliptical distribution 
with density function 



P +iF q ( ai ,...,a p ,a+ I 6 9 ; Y ( S + iS~ 1/2 (X - ^'©-^X - ^S" 1 



/2 



|S+ |S-V2(X - / 1 )'©- 1 (X - /1 )s- i / 2 | a+n/2 

«n£/i constant of proportionality 

r m [a + n/2] |H| a 



(27r)™«/2r m [a]|S|»/2|0|™/2 3 , +1 F a (ai, . . . , a p , a; &i, 
w/iere Re(a) > (m — l)/2. 



,6 9 ;Y E' 1 ) 



Proof. Follows immediately form [3] and Lemma 12.21 □ 
Observe that, by taking Y = in Theorem 1 4.11 we obtain that X h as a matricvar iate 

T distribution, see lDickevI jl967t ). iBox and Tiao 1 (|l972l pp. 441-448) and lPressI (|l982L pp. 

138-141). Also, if we take p = 0, q = 1 and a — bi we obtain that X has a noncentral 

matricvariate T type 2 distribution. Then observing that 



1 F 1 (a;a; Y E~ 



ctr{Y H -1 } 



and by the Kummer relation (jMuirheadl . 1 1 9 8 2| . eq. (6), p. 265) 



1 F 1 [a+ ^;a; Y (s + fS^X - ^'©^(X - ^IT 1 ^ 



= etr \ Y (S + ±5T ^(X - /i)'© _1 (X - A*)S" 1/2 



XiF x a; -Y (s + §5rV2(X - ^'©^(X - m)^ 72 ) * 



Observe that for n/2 an inte ger, iFj is a poly nomial of degree mn/2. In this case the density 
of X is evaluated easily, see iMuirhead (|l982l p. 258). 
From and Lemma \2. 3 1 we have the following result. 

Theorem 4.2. Assume that X|P - Af mXTl (iJ,, S 1 / 2 PS 1 / 2 , ©), where P has a matrix variate 
inverted hypergeometric beta type II distribution. From Lemma \2.3\ its density function is 



</ p (P) oc |P| b -( m+1 )/ 2 |I m +P|-( a+b ) p ^(ai, 



,a p ; HP 1 ), P > 0, 



S 



where the constant of proportionality is 



P m [a,b]p+iF g+1 (ax, ...,a p ,a; bt, . . . , b q , -b + (m + l)/2; — S) ' 

and Rc(a) > (to — l)/2 and Re(6) > (to — l)/2 + fci. T/ien X has a matrix variate left- 
spherical elliptical distribution with density function 



oc p+ iF q+1 (ai,...,a p ,a + n/2; b x , . . . , b q , -b + 



(m + 1) 



- s + ±sr 1/2 (x - mJ's-^x - /x)s- 1/2 

where the constant of proportionality is 

(2tt) ~'"" /2 /3 m [a + n/2, b - n/2] | S | ~"/ 2 1 1 ~ m / 2 

j9 m [a,6]p+iFg + i(ai, . ..,a p ,a+ n/2; 61, . . . , 6 q , -6 + (to + n + l)/2; — S) ' 

where Re(a) > (to - l)/2, Re(b) > (m + n - l)/2 + 

A result of particular interest is obtained from Theorem 14.21 ta king 5 = . Similarly, 
in the Bayesian context, Theorem 14.21 generalises a result given in Etui (jlQQOh . which can 



be obtained by taking 5 = and p = q = 0. In this latter case, by applying the Rum- 
mer relation (Muirhead, Il982l . Theorem 7.4.3, p. 265) we obtain a matricvariate confluent 



hypergeometric of the first kind distribution type. 

Theorem 4.3. Assume that X|P - Afmxniv, £ 1/2 P£ 1/2 , 0), where P has a matrix variate 
inverted generalised hypergeometric distribution. By Lemma[ 



1. its density function is, 

5p (P) cx iPl-^+^ViMsq-SP- 1 ), P > 0, 
where the constant of proportionality is 

l B l"/3 m [a,c-q] 
(3 m [a, b - a]f3 m [a - a, c - a) ' 

with Re(a) > (m— l)/2, Re(a - a) > (m - I)/2, Rc(6 - a) > (m - I)/2, Re(c - a) > 
(to - l)/2 andRe(c-a) > (m - l)/2. 

I?. Or with density function 

, 9p (P) cx |P|-«-("*+ 1 )/ 2 1 f 1 (6;q -SP- 1 ), P > 0, 

where the constant of proportionality is 

\S\ a T m [b}T m [c-a} 
T m [a\T m [c]T m [b - a] ' 

with Re(a) > (to - I)/2, Re(6 - a) > (m - l)/2 and Re(c-a) > (to - l)/2. 

TTien X /ias a matrix variate left- spherical elliptical distribution and 
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1. its density function is 

oc |E _1 (X - m)'© -1 (X - /x)r (a+ ™ /2) 

x 3 Fi fa, 6, a + n/2; c; -2E (lT 1/2 (X - fi)'e _1 (X - /i)£ _1/2 ) 

where the constant of proportionality is 

2 ma \S\ a T m [a + n/2]f3 m [a 1 c-a] 
n mn / 2 (3 m [a,b- a]0 m [a -a,c- a]|E|"/ a |e| m /2 ' 

with Re(a) > (m-l)/2, Re(a-a) > (m-l)/2, Re(b-a) > (m- l)/2, Re(c-a) > 
(rn - l)/2 and Re(c - a) > (m- l)/2. 

I?. Or with density function given by 

oc IXT^X - /x)'0" 1 (X - /x)r (a+ ™ /2) 

x 2 Fi (a,a + n/2;&;-2H(S" 1/2 (X-/2)'0" 1 (X-/x)S" 1/2N 



where the constant of proportionality is 

2 ma \S\ a T m [a + n/2]\S\ a T m [b]T rn [c - a] 



7r"W 2 r m [ a }T m [c}T m [b - a] | S | ™/ 2 1 | m /2 ' 
with Re(a) > (to - l)/2 ; Re(6 - a) > (to - l)/2 and Re(c-a) > (to - l)/2. 
Proof. Follows from ((3]) and Lemma \2. 61 □ 

5 Compound matricvariate T distribution 



From|Dickej (|l967l ). iBox and Tiao I (|l972l . pp. 441-448) and lPressI (|l982l pp. 138-141) we 



know that X has a matricvariate T distribution, denoting this fact as X ~ _A/f7" mxn (z^, /x, X, 0) 
if its density function is 



7r m "/ 2 r m [^/2]|S|»/ 2 |0| m /2 
where 1/ > m — 1. 



Theorem 5.1. Assume that X|P ~ MT mxn (i/, /i, S 1 / 2 ?^ 1 / 2 , 0), where P /ias a matrix 
variate inverted hypergeometric beta type II distribution with 3 = 0. From Lemma \2.3\ its 
density function is 

g p (P) oc \P\ b -( m+1 V 2 \I m +P\-( a+b \ P > 0, 
where the constant of proportionality is 

1 



P m [ a , b Y 

and Re(a) > (m — l)/2 and Re(6) > (m — l)/2. Then X has a matrix variate left- spherical 
elliptical distribution with density function 

„ ( (n + v) n , (m + n + 1) „ 1 1 ,„ 
« 2^1 I L ~^- J ; a +2 ; ~ b+ 2 ' S ~ (X ~ M) ®~ (X ~ M) 
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where the constant of proportionality is 

T m [(n + v)/2](3 m [a + n/2, b - n/2] 
x™/*r m [v/2]P m [a,b]\V\ n ' 2 \e\ m '*' 
where Re(a) > (to - l)/2 and Re(6) > (m + n - l)/2 + ki. 
Proof. Follows from Lemma 12.31 noting that 

-fn+v)/2 

i m + s- 1 /2( X - /xj'e-^x - ^s- 1 /^- 1 

By applying the Euler relation (jMuirheadl . [l982L eq. (7), p. 265) to results in Theorem 15.1 
the density function of X is then 

Vn-l/Y ,|-( a +b+(n+v)/2-(m+l)/2) / i/ m + l 

oc |S (X-/x)0 (X-/x)| 2^1 I -6- 2+^—^ 

m + l , (to + n + 1) „ i, 1, 
- a - 6 + — ^— ; -6 + ^ IT^X - m)'© _1 (X - #») 

where, if a and 6 are integers, 2 (a + 6) > to + 1 and to is odd, 2-F1 is a polynomial of degree 
m(a + 6 — (to+ l)/2). Similarly, if 6 and z//2 are integers, 2 (q + y/2) > m+ l and m is odd, 
2-F1 is a polynomial of degree m(a + v/2 — (to + l)/2), see iMuirhead (|1982L p. 258). 

6 Scale mixture of elliptical distribution of a random matrix 

The a pproach known as the scale mixture of no rmal distributions , prop osed bv iMuirheadl 
( 1982 . p. 33) for the vector case and extended bv lGupta and Varea ( 19931 Chapter 4) to the 
matrix variate case, is obtained as a particular case of the approach described in the Section 
|3l To do so, we take to = 1 in the distribution G(P), from where we ob tain the following 
™ach, termed the scale mixture of an elliptical distribution, cited bv lGur^ and Vargal 
(I 19931 . pp. 78-79). 

Assume that the conditional distribution 

X\s~££S mxn (ii,<t>(s)-E,®,h), (4) 

where <fi '■ (OjOo) — > (0, 00), with s > has the distribution function G(s). Then X has a 
left-spherical elliptical distribution (scale mixture of elliptical distribution) with a density 
function given by 

c(m,n) f ,,_ mn /2, / 1 



f (<K«))~ mn/2 h (^) S_1(X _ ^) ,6rl ( x - dG ^- ( 5 ) 



| S |n/2| |m/2 

As an example, consider the following version of Theorem 12.61 for to = 1. 
Theorem 6.1. Assume that X|s has a matrix variate normal distribution, 

X|s ~ Nrnxnb*, « S 7 ©), 

where s has an inverted hypergeometric gamma type distribution. By Lemma \2.S\ its density 
function is 



g s (s) ocexpj-^js a l p F q (a,!, . . . , a p ; b x , . . . , b q ; ^ 



s > 0, 
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where v > 0, £ > and the constant of proportionality is 



and Re(a) > 0. TTien X has a matrix variate left- spherical elliptical distribution with density 
function 

p+ iFq (ai,. ..,a p ,a+ A, ■■■,b g ;v(£ + \ trE _1 (X - /x)'0 _i (X - fj,)) ^ 
^ (£ + I trS-i(X - M )'0-i(X - M )) a +™/ 2 

constant of proportionality 

r[a + ron/2] £ a 



(27r)™'V2r[a]|S|'V2|0|-/2 p+lj p 9 (ai, a p , a; &i, ~) 
where Re(a) > 0. 

Similarly, from Theorem 16.11 particular cases are obtained, taking, for example, v = 0, 
but in this case we obtaining the matrix variate T distribution (not the matricvariate T 
distribution). Similar consequences are obtained as particular cases, taking m = 1 from the 
Theorems S2K5II1 

Conclusions 

This paper introduces several hypergeometric type distributions, which include many that 
are well-known in the statistical literature, such as central and noncentral matrix variate 
inverted Gamma (Wishart) and matrix variate beta type II distributions, among many 
others. 

Assuming that P has one of these hypergeometric distributions in 
X|P ~ ££S mxn (n, S 1 / 2 *(P)S 1 / 2 , 0, h), 

several left-spherical matrix variate elliptical families are found. These enable us to study 
examples in which the tails of the distributions are heavier or lighter than in the normal 
case. These approaches to obtaining elliptical distributions ar e of particular interest from the 
Bay esian standpoint and for s hape theory distributions, see Ijammalamadaka et all (|l987t l 
and ICaro-Lopera et al. ( 2008 ). respectively. 
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